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Abstract. 

. . . We present the first results of a study where we determine the metallicity distribution function 

OO in the Galactic disks as a function of height above the Galactic plane. Observations in the 

^^ Stromgren photometric system enables us to identify the dwarf stars and derive metallicities 

for them. The resulting metallicity distribution functions at 0.5 and 2.0 kpc above the Galactic 
plane are significantly broader and more metal-rich than is anticipated from standard models 
such as the Besangon model. Our results can be explained by invoking a smaller scale height 
and larger local normalisation for the thick disk than is commonly used in the models. These 
results are compatible with recent determinations of the thick disk scale height based e.g. on 
SDSS data. 

The age of the stellar populations as a function of height above the Galactic plane is also 
investigated by studying the turn-off colour and metallicity. We tentatively find that at 2.0 kpc 
,^^ 1 above the Galactic plane there exist an intermediate age population. 

1 Keywords. Galaxy; stellar content - Galaxy: structure - Galaxy: disk - stars: distances 

O 

C^ . 1. Introduction 

The Milky Way galaxy hosts two stellar disks. The properties of the disks appear to 
differ in that the thick disk is older and more metal-poor than the thin disk. Much of our 
r\ ' knowledge of these properties for the disks is acquired through the study of local stellar 
^ ; samples (e.g. Bensby et al. 2005, Fuhrmann 2008, Reddy et al. 2006). 
\^ The stellar sample studied here lies along the line-of-sight towards the Draco dwarf 

T— I spheroidal galaxy, located at I — 86.4 and b = 34.7, Fig. [1] Along this line-of-sight the 

r — , Milky Way disk stars are approximately at the same galactocentric distance as the sun. 
f^ ■ This means that we are fairly insensitive to any radial gradients that might exist in 
OO , the disk, and can assume that any gradients found are basically due to changes in the 
^^ properties of the stellar populations in the z-direction. 

^ , This study uses the Stromgren photometric system which is made up of four intermedi- 

l^ ' ate band filters (uvby) that are strategically located so that information can be obtained 
i*^ , about the amount of line blanketing and the size of the Balmer jump. Stromgren (1963) 
C^ ■ devised two indices, mi and ci, that measures these effects. While the mi index is metal- 
licity sensitive, the ci index is sensitive to surface gravity and we use it to separate 
distant giant stars, belonging to the dwarf spheroidal galaxy, from Galactic dwarf stars 
and sub-giant stars, belonging to the thin disk, thick disk, and halo, respectively. 

2. The data and parameters derived from the photometry 

The observations were obtained with the Isaac Newton Telescope on LaPalma using 
the Wide Field Camera. Observations were taken both of the centre and the outskirts of 
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Figure 1. Schematic figure showing the observed hne-of-sight (LOS). 

the Draco dwarf spheroidal galaxy. After carefully separating the background members 
of the dwarf spheroidal galaxy from the foreground dwarf stars and sub-giant stars we 
present here results based on the properties of 749 Galactic dwarf and sub-giant stars 
with y == 15 - 18.5. 

For these stars we have determined both distances and metallicities based on Stromgren 
photometry. The metallicity calibration of Ramirez & Melendez (2005) was used for 
stars with {b — y) ^ 0.8 and the metallicity calibration of Olsen (1984) for stars with 
0.8 < {b — y) ^ 1.0. Of the available calibrations in the literature these best reproduce 
spectroscopic iron abundances for a compiled set of 460 test stars (Arnadottir et al. , in 
prep). In this magnitude range the errors in the derived metallicities are below ^ 0.3 
dex. 

The line-of-sight toward the Draco dwarf spheroidal galaxy has been observed as part 
of the Sloan Digital Sky Survey (SDSS) so ugriz photometry is available for our sample 
of stars. We have estimated the completeness of our photometry using the SDSS survey, 
and base our magnitude cut on that completeness estimate. The distance modulii of indi- 
vidual stars were determined using ugriz apparent magnitudes and absolute magnitudes 
estimated using the isochrones of Girardi et al. (2004). For stars near the turnoff we 
calculated a logg using our own calibration (Arnadottir et al., in prep). This logg was 
used to constrain the solution, i.e. we know if we were determining the distance to a main 
sequence dwarf or a sub-giant star. 

The photometric ugriz metallicity calibration of Karaali et al. (2005) was tested but it 
clearly is not able to derive the metallicities for the dwarf stars under consideration here 
and can therefore not be used for a larger investigation of the dwarf star properties in the 
disk based on ugriz photometry. Additionally, there is very little overlap between this 
study and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) 
survey. Only 4 stars were found in common between the two. 



3. Metallicity distribution as a function of height above the Galactic 
plane and implications for the scale heights of the thin and the 
thick disk 

We present the metallicity distribution functions at 0.5 kpc and 2.0 kpc above the 
Galactic plane, Fig.[2l At both heights we see a broad distribution of metallicities. 

The observed distributions are not well reproduced by the Besangon model of the 
Galaxy (Robin et al. 2003). The predicted metallicity distribution from the Besangon 
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Figure 2. The metallicity distribution function for dwarf stars with 15 < 1/ < 18.5 is shown at 
heights a) z=0.5 kpc (±0.2 kpc) and b) z=2.0 kpc (±0.5 kpc) above the Galactic plane in the 
direction of {I, h) = (86.4, 34.7). Typical error in [Fe/H] is < 0.3 dex. 



model is narrower and at 0.5 kpc above the Galactic plane it is more peaked towards 
solar metallicities and the model predicts that 80% of the stars are contributed from the 
thin disk and 20% of the stars are from the thick disk. At 2.0 kpc above the Galactic 
plane the predicted metallicity distribution from the Besangon model is also narrower 
than the observed one, and peaks at lower metallicities, [Fe/H]~ —0.9. Here the model 
predicts that 89% of the stars are thick disk stars with only 7% contribution from the 
thin disk and 4% from the halo. 

In summary, our data show broader metallicity distribution functions than the Be- 
sangon model and our data also peaks at higher metallicities than the model does. Can 
we reconcile our data with a realistic model of the Milky Way disk system? 

One thing that is not too well constrained in the literature are the scale heights of the 
thin and the thick disks and the local normalisation of the stellar population densities. 
The thin disk scale height is determined to be between 220 pc and 330 pc (see i.e. Du et 
al. (2006) and Chen et al. (2001)). The thick disk scale height, varying between 580 pc 
and 1600 pc, and its local normalisation, varying between 1% and 13%, are even more 
unconstrained. In Fig. [3] we have compiled a collection of thick disk scale heights and 
local normalisations from the literature. There is a large range of acceptable values. The 
BesauQon model (Robin et al. 2003) has a thick disk scale height of 800 pc and a local 
normalisation of 6.8%. 

Using a toy model where each stellar population, i.e. thin disk, thick disk and halo, is 
represented by a single age and single metallicity we varied the scale heights of the thin 
and the thick disk and the local normalisation. This enabled us to study the contribution 
from stars (15 < V < 18.5) of each population to the metallicity distribution at 0.5 kpc 
and 2.0 kpc above the plane, and qualitatively compare it to the observed metallicity 
distributions shown in Fig. [51 We tested several combinations of thin and thick disk 
scale heights and local normalisation from the literature (Fig. [3]). The observed broad 
metallicity distributions, Fig.[2l suggest that the contribution from the thin and the thick 
disk are much closer to 50-50% than the model predicts. At this stage our comparison 
can only be qualitative. However, we believe that the results will indeed hold up also 
with the more detailed modelling that we are now undertaking. 

Our preliminary results from this simple modelling is that a thin disk scale height of 



Anna S. Arnadottir, Sofia Feltzing & Ingemar Lundstroni 



O Gilmore & Reid (1983) 

D Gilmore & Wyse (1985) 

A Robin & Creze (1986) 

D Kuijken & Gilmore (1989) 

A Reid & Majewski (1993) 

V Robin el al. (1996) 
O Spagna et al. (1996) 

V Norris (1996) 

X Ng el al. (1997) 



Buser el al. (1998) 

Buser el al. (l999) 

Chen al. (2001) 

Kerber el al. (2001) 

Ojha (2001) 

Reyle & Robin (2001) 

Du el al. (2003) 

Larsen & Humphreys (2003) 

Du el al. (2006) 



0.15 


- 


1 




' 




' 


1 




1 


' 




ization 

o 




•. 


















; 


Id 

E 

u 
o 

c 

13 0.05 




-^■♦. 


'(■ 

■■■■.V 


T 


-)f 


'^:::: 


■::,.,..+. 


■v 






- 


o 

q 


- 








A-.. 


D 








- 


n 






1 1 




-k 


"d-.-.- 


^■■x ^"■"■"- 




.■..^.vQ.-V.v. 


■ ./A 





600 



800 1000 1200 

Thick disk scale height 



1400 



1600 



Figure 3. The local normalisation is plotted against the thick disk scale height for a compiled 
list of values from the literature. References are given above the Figure. When a study gives a 
range of values, the extremes are plotted and connected with a dotted line. 

320 — 320 pc, a thick disk scale height of 580 — 640 pc, and a local normalisation of 
7 — 13% appear to explain our data. This combination of scale heights is advocated by 
recent studies such as Chen et al. (2001), which is based on star counts in 279 deg^ in 
the SDSS survey, and Du et al. (2006), which are based on star counts in 21 fields of 1 
deg^ in the Beijing- Arizona- Taiwan-Connecticut (BATC) survey. 

In order to try to assess whether the observed line-of-sight is representative for the 
Milky Way disk system we queried the SDSS for stars in eleven fields in the northern 
galactic hemisphere with 90 < Z < 270 and 35 < 5 < 90 (including one field centred 
on the Draco dwarf spheroidal galaxy). We systematically compared colour magnitude 
diagrams and the colour distributions in these fields to the results from the Besangon 
model. 

We find that in each field and in the magnitude range 15 < V < 18.5 the turn-off stars 
in the Besangon model are systematically bluer and overproduced as compared to the 
SDSS data. In this magnitude range we are, in the bluest part of the colour distribution, 
looking at the thick disk turn-off stars at ~ 1 — 3 kpc above the Galactic plane. The 
observed overproduction is in accordance with our previous indications that the true 
scale height of the thick disk is thinner than that used by the Besangon model. 

Since this effect is observed systematically along several lines-of-sight we conclude that 
the line-of-sight towards the Draco dwarf spheroidal galaxy is not unique and conclusions 
drawn here about the thick and thin disk scale heights from the metallicity distributions 
along that line-of-sight are of a more general nature and can be applied to the Milky 
Way as a whole. 



4. Population age as a function of height above the Galactic plane 

We are able to study the turnoff colour as a function of height above the Galactic plane, 
enabling us to measure the age of the population and their metallicities as a function of 
height. In Fig. [4] we show the colour distribution of dwarf stars at different heights above 
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Figure 4. Metallicity is plotted against ib — y) colour for stars with e[Fe/ir] < 0.5 and log <; > 3.8 
at three heights above the Galactic plane: a) 1.0 ± 0.3 kpc, b) 1.5 ± 0.3 kpc, and c) 2.0 ± 0.5 
kpc. Dashed lines show turn-off colour and metallicity for three ages (2 Gyrs, 8 Gyrs and 14 
Gyrs). The turn-off colours of the globular clusters M92, NGC288 and 47 Tuc are marked with 
diamonds. 

the Galactic plane. Only dwarf stars are plotted. The turnoff colour and metallicity are 
shown for 3 sets of ages, obtained using the isochrones of VandenBerg et al. (2006). By 
determining where the turn-off stars lie, one can read off the population age at that 
height above the Galactic plane. For validation we mark three globular clusters, M92 
([Fe/H]= -2.14), NGC288 ([Fe/H]= -1.31), and 47 Tuc ([Fe/H]= -0.70) that can aU 
assumed to be old, onto the plot (Calamida et al. 2007). 

Fig. [H shows an indication of an intermediate age population (younger than 8 Gyrs) 
present even at large heights above the galactic Disk in the direction towards the Draco 
dwarf spheroidal galaxy (the direction our observations were taken in) . 



5. Summary 

Using Stromgren photometry to derive the stellar parameters of dwarf stars in the 
Milky Way along the line-of-sight towards the Draco dwarf spheroidal galaxy we have 
studied the properties of the stellar populations in the Milky Way disks and find that 

• The metallicity distribution functions at 0.5 kpc and 2.0 kpc above the Galactic 
plane are broader than predicted by the Besangon model and the observed metallicity 
gradient as a function of height above the Galactic plane is not as steep as predicted. 

• Preliminary investigations, using a simple toy model, indicates that our observations 
favour models with smaller thick disk scale heights and higher local normalisation such 
as found by e.g. Chen et al. (2001). 

• By studying the turn-off colour and metallicity as a function of height above the 
Galactic plane, we find an indication of an intermediate age population located 2 kpc 
above the Galactic plane. 

Further modelling of the data is being undertaken, but we believe that these prelimi- 
nary results will hold. 
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